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ABSTRACT 



CO I Luminosity functions and their integrated luminosity densities are presented for 

the 6dF Galaxy Survey (6dFGS). This ongoing survey ultimately aims to measure 
around 150000 redshifts and 15 000 peculiar velocities over almost the entire southern 
sky at |6| > 10°. The main target samples are taken from the 2MASS Extended Source 
Catalog and the SuperCOSMOS Sky Survey catalogue, and comprise 138 226 galaxies 
^ ; complete to {K ,H ,J ,r^,b,) = (12.75, 13.00, 13.75, 15.60, 16.75). These samples are 

\^ , comparable in size to the optically-selected Sloan Digital Sky Survey and 2dF Galaxy 

ff^ ■ Redshift Survey samples, and improve on recent near-infrared-selected redshift surveys 

I by more than an order of magnitude in both number and sky coverage. The partial 

. samples used in this paper contain a little over half of the total sample in each band 

' and are ~90 percent complete. 

fH \ Luminosity distributions are derived using the 1/V^nax, STY and SWML estima- 

tors, and probe 1 to 2 absolute magnitudes fainter in the near-infrared than previous 
surveys. The effects of magnitude errors, redshift incompleteness and peculiar veloc- 
ities have been taken into account and corrected throughout. Generally, the 6dFGS 
luminosity functions are in excellent agreement with those of similarly-sized surveys. 
Our data are of sufficient quality to demonstrate that a Schechter function is not 
an ideal fit to the true luminosity distribution, due to its inability to simultaneously 
match the faint end slope and rapid bright end decline. Integrated luminosity den- 
sities from the 6dFGS are consistent with an old stellar population and moderately 
declining star formation rate. 



o 



1 



Key words: surveys — galaxies: clustering — galaxies: distances and redshifts 
cosmology: observations — cosmology: large scale structure of universe 



1 INTRODUCTION 

Arguably the most fundamental of all cosmological observ- 
ables is the mean space density of galaxies per unit lu- 
minosity, or luminosity function (LF). In the local uni- 
verse, the shape of the LF places strong constraints 
on models of dark halo formation and galaxy evolu- 
tion. Current semianalytic models that invoke feedback 
in the form of supernovae and stellar winds are un- 
able to satisfactorily reproduce the sharp bright-end de- 
cline, while hydrodynamic simulations typical l y produce 
LFs t hat are too bright ('e.g. iBenson et al]|2003l: iKav et all 
|2002| . and references therein). Alternative approaches (e.g. 



Coorav fc Milosavlievid 120051: iBerUnd fc Weinberg l2002l : 
^ng^M^^^ar^Bir^oschll200a) use the LF as an empir- 
ical constraint when inferring how galaxies populate dark 
matter halos, thereby providing testable predictions about 
galaxy clustering. 

In the near- infrared (NIR), LFs are effective tracers of 
the stellar mass function for many reasons. Unlike the visi- 
ble or far-infrared, light from the NIR is dominated by the 
older and cooler stars that make up the bulk of the stellar 
mass. Compared to a survey at visible wavelengths, the bal- 
ance of galaxy types in a NIR survey shifts from late to early 
types, which contain most of the stellar mass of the universe. 
Furthermore, the effects of extinction are minimal at longer 
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wavelengths, meaning that NIR luminosities are affected lit- 
tle by dust, either in our own Galaxy or in the target galaxy. 
Mass-to-light ra tios are also much be tter constrained in the 
NIR passbands llBell fc de Jond lioOll. 



Our current understanding of the galaxy LF owes 
much to the 2dF Galaxy Redshi f t Survey (2dFGRS; 



Folkes et alJ Il99(5l: ICross et alJ I2OOII: iNorberff et alJ 120021 


Mads 


rwick et al.l I2OO2I: 


De Prooris et al.l 120031: lEke et alJ 


20o4 


ICroton et alJl2005 


) and the Sloan Digital Sky Survey 


(SDSS;"Blanton ct al."200l'; 'Goto ct al. 2002; 


Blantoi^^l 


2003; 
I2OO5I1 


Baldrv et al., 2005 


Berlind et al. 2005: 


Blanton et al. 


. The wide sky coverage and extremely large optically- 



selected samples (~ 10^) utilised in these studies have 
produced a tenfold increase in the precision with which 
the shape of the LF is known. With t he advent of 
the Two-Micron A ll-Sky Survey (2MASS: I.Tarrett et alJ 
| 2000l). other studies iCole et aLlbOoltlKochanek et alJl200lf 
iBell et al ."2003':'E ke et alJl2005^ have combined the 2MASS 
NIR galax y samples with the 2dFGRS, SDSS and ZCAT 
iHuchra et al.lllS?93) redshift catalogues to yield wide sky 
coverage and moderately large NIR-selected samples ('^ 
10^). 

In this paper we present KHJrYb,, LFs from an interim 
subse t of the 6dF Galaxy Survey (6dFGS: Jone s ct al. 2004,, 
|2Qq3). a NIR-selected redshift survey of 150 000 galaxies 
over almost the entire southern sky combined with a 15 000- 
galaxy peculiar velocity survey. The full target list contains 
174 442 sources, although most of these (138 226) come from 
five main samples that are complete to {K , H , J ,rF,bj) = 
(12.75,13.00,13.75,15.60,16.75). The remainder are pro- 
vided by a do zen supporting samples selected in various ways 
l|.Tones et al.l Enn4'l. When complete, the 6dFGS will cover a 
sky area eight times that of 2dFGRS and twice that of the 
revised SDSS. The median 6dFGS redshift (2 = 0.054) is 
roughly half of that of the other two surveys, and reflects its 
shallower limiting magnitudes. At completion, the 6dFGS 
will be two-third s the sample size of 2dF GRS and about 
one-fifth of SDSS. I Jones et al.l ||200l|200i) give full descrip- 
tions of the survey and show redshift maps of the galaxy 
distributions obtained thus far. 

As well as providing initial results on the 6dFGS LFs, 
this paper gives a comprehensive explanation of the LF 
derivations used in this and related future papers. Section 2 
summarises the target selection, galaxy photometry and cor- 
rections for incompleteness; surface brightness selection is- 
sues are also discussed. Section 3 examines normalisation 
issues via number count data, and explains how we have 
accounted for galaxy magnitude errors, peculiar velocities, 
extinctions and fc-corrections. In Section 4 we present the 
6dFGS LFs and compare the results to other recent work. 
Section 5 is devoted to a comparison of the 6dFGS luminos- 
ity density across b.^rpJHK with other relevant surveys. 

We have adopted the current standard cosmology 
(JIm ~ 0.3, f^A ~ 0.7) throughout the paper. 



2 TARGET SELECTION AND PHOTOMETRY 

2.1 Areal and Magnitude Selection 

The 6dF Galaxy Survey commenced observations in 2001 
May although formal observing did not start until 2002 Jan- 



uary. Incremental public data releases were made in 2002 De- 
cember, 2004 March, and 2005 May, the latter making avail- 
able on-line 83 014 sources with their spectra, redshifts and 
photometric measurements. Ijones et alJ ll2004) describe the 
First Data Release (2004 M arch) and the char acteristics of 
the 6dFGS in general, while Ijones et all ll2005l) describe the 
updates in place for the Second Data Release (2005 May). 

In total, there are 1564 fields and 174 442 galaxies in 
the 6dFGS, yielding a mean source density of 112 galaxies 
per field. With an instrument field size of 25.5 deg^ and 
10.2 sources deg~^, the field tiling has effectively had to 
cover the southern sky twice over. FigureQshows the survey 
coverage to date from the 1 390 fields contributing to the LFs 
herein. Although substantial areas at mid- and equatorial 
declinations have been covered, some areas await a second 
pass, and as a result, redshift completeness will ultimately be 
higher. Coverage around the Pole is still largely incomplete, 
but will be addressed during the closing stages of the survey. 

Targets for the 6dFGS were selecte d from the 2MASS 
Extended Source Catalog (2MASS XSC: IJarrett et al.l200o)) 
for southern declinations 5 < 0° and at galactic latitudes 
&| > 10° for JHK and |6| > 20° for 6.,r-F. Table □ sum- 
marises the selection criteria, sample sizes and sky coverage. 
T he motivations for t he sample selection criteria are outlined 
in I Jones et al] Jiooi). 

Total magnitudes are the preferred way of quantifying 
galaxy fiux because the related observable (total luminosity) 
is more physically meaningful than isophotal or aperture 
magnitudes that depend on selection criteria or distance. 
At the time the target lists were prepared, however, the 
total K magnitudes (-f^tot) from the 2MASS XSC were not 
sufficiently reliable at the lowest galactic latitudes in the 
survey due to insufficient depth and spatial resolution. The 
isophotal 2MASS K magnitudes (A"iso) were more robust, 
however, and so formed the basis for a derived Ktot that we 
used for target selection in K: 

Ktot = ifiso - 1.5 exp [1.25{TrW20 - 20)]. (1) 

Here, fiK20 is the mean surface brightness within the hk = 
20 elliptical isophote. The scatter in (-K^iso — ^tot) is typically 
no more than 0.1 to 0.2 magnitudes over the r ange of surface- 
brightnesses encountered (Fig. 1 of l,Tones~ al. 2004 '). 

This correction was eventually superseded by the re- 
vised total K magnitudes from 2MASS XSC which became 
available after 6dFGS observations had commenced. These 
show less scatter than the Ktot values used for the 6dFGS 
sample selection (T. Jarrett, priv. comm.) and are a superior 
measure of total flux. 

Both the old and new magnitudes contribute to the 
6dFGS luminosity functions in different ways. The origi- 
nal Ktot magnitudes are preserved for selection purposes 
and were used for deriving the selection function and com- 
pleteness corrections (Sect. 12.311 . However, the superior K 
magnitudes are used to compute galaxy luminosities after 
all selection cuts and completeness corrections have been 
made. 

The first studies to use 2MASS pho tometry for LF 
work were divided over the use o f isophotal llKochanek et alJ 
[ 200 ll) or total magnitudes JCole et alJ I2001II . While 
iKochanek et alJ claimed their (j ^iao — -ft'tot) correction was 
only 0.20 mag, lAndreoiil l|2002) argued that it is closer 
to 0.3 mag. Claims that differences between the optical 
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Figure 1. (a) Equal area Aitoff projection of the total 6dFGS field coverage (circles) superposed with the subset of fields contributing 
to this paper (filled circles), (b) Galaxy density (per square degree) of 6dFGS targets on the sky. 



SDSS and 2MASS NIR luminosity densities were d ue to a 
2MA SS bias against low surface brightness systems iWrighli 
120011) were largely countered when the SDSS optical lu- 
minosities were re-c omputed to include galaxy evolution 
jBlanton et alJl2003l) . although a small offset remains. 

ICole et all ll200ll) have undertaken a comparison of 
their adopted 2MASS Kron and extrapol ated magnitude s 
with the deeper infrared photometry of iLovedavl i200Cl) . 
They found scatter of around cr ~ 0.13-0.14 magnitudes in 
both samples, although found that this reduced to a « 0.1 
if K Kron magnitudes were derive d from higher qu ality J 
data and (J — K) colours. Likewise. iBell et al.l i2003l) find a 
scatter if ~ 0.2 mag in comparisons of 2MASS Kron magni- 
tudes with those of Lovcdav (200^. Since these studies, the 
2MASS photometry has been revised. With the current to- 
tal magnitudes, we find (i^iso — -R'tot) = 0.08 for A'iso < 10, 
rising to (Aiso — Atot) = 0.12. At A'iso 11, the correc- 
tions can be as large as 0.5 mag, reflecting the lower sur- 
face brigh t ness l ate-type galaxies that populate this regime. 
lEke et al.1 ^2005^ find that the improved 2MASS total mag- 
nitudes match the lLovedavl data with a scatter of 0.125 mag. 

The original bjrp samples for 6dFGS were selected from 
the SuperCOSMOS catalogue in 2001. This photometry suf- 
fered from scatter in the fleld-to-field zeropoints of a few 
tenths of a magnitude. After the all-sky digital photometry 
of 2MASS subsequently became available, the original 6,,, rp 
(and jn) SuperCOSMOS magnitudes were re-calibrated in 
2003 to a common zeropoint by J. Peacock, N. Hambly and 
M. Read, for the benefit of the 2dF Galaxy Redshift Sur- 



vey. Individual field zeropoints were characterised by the 
mean (6,i — J) colour. The revised SuperCOSMOS mag- 
nitudes show a significant improvement in the photometric 
calibr ation, with zerop o int scatter of 0.04 mag rms in each 
band iCole et alJl2005h . ICross et al.l J2004ft find a scatter of 
(T = 0.1 mag in a comparison of the revised SuperCOSMOS 
magnitudes with the CCD photometry of the Millennium 
Galaxy Catalogue and Sloan Early and First D ata Releases 
JStouehton et al.ll2002l: I Abazaiian et al.ll2003D . 

As with the revised K photometry, the original fejrp 
magnitudes were retained for selection and completeness 
calculations while the re-calibrated ones were used for lumi- 
nosities. This is similar to the approach used by the 2dFGRS 
when re-calibrated 6j photometry became available mid-wa y 
through that survey jNorberg et alj|2002l: fCole et al.ll2005l) . 
However, in that case, the new and old magnitudes were 
related by a simple transformation which imparted no ad- 
ditional scatter. In the case of 6dFGS, the relationship be- 
tween the old and new magnitudes is not as straightforward. 
However, the scatter that results is small and inconsequen- 
tial, largely because the dataset itself remains the same. 



2.2 Surface Brightness Issues 

One common misgiving surrounding apparent magnitude 
limited surveys is the extent to which t hey bias the 
LF against low surfac e brightness systems iPisne ^^T^: 
Imp ev fc BothunI Il997l: ISoravberrv et alJ Il997: Da lcantonI 
[1993). Analyses utilising bivariate brightness distributions 
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Table 1. 6dFGS target samples 



Sample Selection 




Full Sample 


This Paper 




Total 


Additional Area 


Number Effective Area 




Number 


to K sample (deg^) 


Used (deg-^) 



2MASS K < 12.75 


113 988 




17045 


60 689 


9 075 


2MASS H < 13.00 


90 407 


3 282 


17 045 


53 142 


10 019 


2MASS J < 13.75 


93 925 


2 008 


17 045 


55 236 


10 024 


SuperCOSMOS rp < 15.60 


66 904 


9199 


13 572 


37156 


7 537 


SuporCOSMOS bj < 16.75 


64138 


9 749 


13 572 


37128 


7 857 



find < 1 mag arcsec"^ change in the peak surface bright- 
ness for the bj-selecte d 2dF Galaxy Redshift Survey across 
observed luminosities ross et al .i200ll:lDriveJ2004l) . while 
ICross fclDrivc r (2002) claim that the faint-end slope is ro- 
bust to surface brightness effects. 

In the case of the 2MASS XSC. ICole e t al.' ("20011) have 
modell ed the select i on ch aract eristics follo wing t he meth- 
ods o f ICross et al.l l|200li and ' Phillipps~'Dayies fc DisnevI 
ilQQOl) . and find no selection bias against low surface bright- 
ness galaxies by 2MASS. Specifically, they demonstrate that 
their galaxies lie well within the theoretical maximum vol- 
umes corresponding to their apparent magnitude and sur- 
face brightness limits (J = 14.45 and 20.5 mag arcsec"^ 
respectively). The claims bv lBell et al.l ll2003ll that surface 
brightness selection impacts the K-hand LF are based on 
brighter surface brightness limits (/ix — 17 to 20) than the 
fiK ^20 limit of the 6dFGS. 

Taken together, these results imply that at the mag- 
nitudes used by the 6dFGS, 2MASS is virtually unaffected 
by bias against low surface brightness galaxies . The 6dFGS 
magnitude limits are ~ 1.5 mag brighter than the incom- 
pleteness regime for the 2MASS XSC, ensuring an effective 
compromise between sensitivity and robust selection for the 
redshift survey sample. 



2.3 Selection Function 

The probability of obtaining a redshift for a galaxy is ex- 
pected to depend on its apparent magnitude K: the fainter 
it is, the less the chance that its redshift will be secured. 
However, in a survey such as 6dFGS, the basic observable 
units are not the galaxies, but the Selds to which they be- 
long. Therefore, the probability that a galaxy redshift will 
be obtained depends not only on its apparent magnitude, 
but also on the overall redshift success rate of its field. 

The survey selection function gives the survey redshift 
completeness as function of both magnitude and sky posi- 
tion. In short, it is obtained by taking each 6dFGS field and 
fitting a functional form to the decline in completeness with 
magnitude (CoUess et al, .2001,; Norbcrg ct al. 2002). This 
form is then normalised to reproduce the overall redshift 
completeness at that point in the sky. 

Before we can compute the survey selection function we 
must first determine the general form for completeness as a 
function of magnitude for each band KHJrpb,,, modulo field 
completeness. Figure |5| shows how this varies if we divide 
galaxies into one of four groups: 
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completeness 
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B mognitude 

Figure 2. Generalised field completeness and its dependence on 
apparent magnitude (solid points with error bars) for fields in the 
ranges (a) 0.9 < Rf ^ 1-0, (set F = 1), (b) 0.8 < Rf ^ 0.9, {F = 
2), (c) 0.7 < Rf ^ 0.8, (F = 3), and (d) 0.6 < Rf ^ 0.7, (F = 4). 
The magnitude bin widths (dotted histogram) and individual fit 
components (dashed curves) are also shown. 



(i) galaxies from high quality fields with completeness in 
the range 0.9 < Rf ^ 1.0 (denoted set F = 1); 

(ii) galaxies from lower quality fields with 0.8 < Rf ^ 0.9 
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(iii) galaxies from mediocre fields with 0.7 < Rf ^ 0.8 
(^ = 3); 

(iv) galaxies from poor fields with 0.6 < Rf ^ 0.7 (F = 

4). 



R{f) 



The completeness, R{f), of field / was determined as 



NM) 



Nutif) + NedFif) 



Nutif) + iV6dF(/) + Ntif) + N,{f) ■ 



(2) 



Here, N^(f) is the number of extragalactic redshifts ob- 
tained in the field, /, and Np{f) is the number of confirmed- 
extragalactic and yet-to-be-observed sources from the par- 
ent catalogue. N{{f) and Nr{f) are the number of redshift 
failures and sources remaining to be observed, respectively. 
Nnt if) and A'edP (/) count the literature and 6dF- measured 
redshifts within each field boundary, respectively. Around 
one-eighth of our redshifts come from the liter ature, mainly 
Huchr a's ZCAT Catalog llHuchra et alJll992^ or the 2dF- 
GRS JColless et al.ll2001^ in the ratio 4 : 1 respectively. 
While observing for 2dFGRS is complete, ZCAT continues 
to be suppleme nted by the Two M icron All-Sky Redshift 
Survey f2MRS: iHuchra et alj|2006ll . which to date has ob- 
tained around 20 k redshifts to is' = 11.25 across the full 
sky. Sources for which spectra were obtained, and which 
were found to be stars, planetary nebulae, Hll regions, or 
other Galactic objects, were excluded. 

At b, < 16.0 {K < 12), ZCAT redshifts represent more 
than one-third of the sample. As a consequence, complete- 
ness below the 2MRS limit is comparatively high while the 
6dFGS continues to fill out completeness at fainter magni- 
tudes {K > 11.25). This dual redshift contribution is evident 
in our measurements of completeness as a function of mag- 
nitude, which exhibit a bump in the declining portion of the 
curve (Fig.|5J. 

The generalised completeness curves were fit by a 
double-exponential, 

Cpijin) — /3f max[0, 1 — exp (m — ^f)] 

-f(l-/3F)max[0, l-exp(m-^F)]. (3) 

The motivation for this m odel was the single exponential 
form adopted by 2dFGRS ijCoUess et alJ boOlf). modified to 
a double-exponential to fit the shoulder. Table |5| lists the 
fit parameters /3f, Mf and according to passband and 
field group {F = 1, 2, 3 or 4). The shoulder occurs at ^f val- 
ues typically 0.5 to 0.75 mag brighter than the survey limits. 
In the near-infrared bands, the initial decline in complete- 
ness is steeper than Eqn. can match. We experimented 
with fitting a single exponential model and in Sect. I4.l1 we 
demonstrate that the luminosity functions are unaffected by 
deviations of this size from the fit. Individual uncertainties 
in Fig.|5|were computed as follows. Given that completeness 
is the ratio of Nz successful redshifts to sources measured, 
then the uncertainty in completeness per bin is 



AC: 



y^N{Nz+2){N~Nz + 1) 
N{N + 3) 



(4) 



With the form of the magnitude-dependent complete- 
ness known, we are in a position to determine total com- 
pleteness, TF{0,rn), as a function of both sky position 



Table 2. Magnitude-dependent completeness parameters Pa- 
rameters describing the best analytic fit to the magnitude- 
dependent completeness in Eqn. Ol . 



Field Completeness Field /3f /^f 









Set F 






(mags) 


(mags) 


X-band: 
















0.9 < Rf 




1.0 


1 


0, 


.80 


14.84 


12.03 


0.8 < Rf 




0.9 


2 


0, 


.74 


14.72 


11.96 


0.7 < Rf 




0.8 


3 


0, 


.71 


14.90 


12.03 


0.6 < Rf 




0.7 


4 


0, 


.70 


15.02 


12.20 


H-band: 
















0.9 < Rf 




1.0 


1 


0, 


.85 


14.66 


12.27 


0.8 < Rf 




0.9 


2 


0, 


.79 


14.60 


12.20 


0.7 < Rf 




0.8 


3 


0, 


.70 


15.12 


12.39 


0.6 <Rf 




0.7 


4 


0, 


.87 


14.24 


12.37 


J-band: 
















0.9 <Rf 




1.0 


1 


0, 


.83 


15.56 


12.99 


0.8 <Rf 




0.9 


2 


0, 


.78 


15.32 


12.97 


0.7 < Rf 




0.8 


3 


0, 


.79 


15.14 


13.11 


0.6 < Rf 




0.7 


4 


0, 


.84 


15.20 


12.87 


rp-band: 














0.9 <Rf 




1.0 


1 


0, 


.92 


16.96 


15.04 


0.8 < Rf 




0.9 


2 


0, 


.86 


16.96 


14.92 


0.7 < Rf 




0.8 


3 


0, 


.83 


17.25 


15.16 


0.6 < Rf 




0.7 


4 


0, 


.71 


20.00 


14.92 


bj-band: 
















0.9 < Rf 




1.0 


1 


0, 


.87 


18.29 


16.10 


0.8 <Rf 




0.9 


2 


0, 


.81 


18.47 


16.22 


0.7 < Rf 




0.8 


3 


0, 


.78 


18.54 


16.22 


0.6 < Rf 




0.7 


4 


0, 


.70 


21.29 


16.40 



and magnitude m. We assume it is separable function with 
the general form 



TF{e,m) = s{e) ■ CfM 



(5) 



for all galaxies in a given field group F. S{9) is a constant 
scaling the completeness of the individual field to that of the 
total completeness in the same part of sky. 



S{6) = R{e) 



■''"bright ^ ' 



dm 



-''"bright ^ ' 



N{m) dm 



(6) 



Actual sky completeness, R{0), is the ratio of extra- 
galactic redshifts {cz ^ 750kms~^) at sky position 0, to the 
number of confirmed or potentially extragalactic sources at 
the same position. It is measured in exactly the same man- 
ner as field completeness R{f), except on a uniform 1.5° 
grid of circular regions of diameter 2°. This geometry was 
unrelated to the placement of actual 6dFGS fields or the 
way in which they overlapped. Figure |3 shows the redshift 
completeness of the all KFUryb.^ samples. Typically there 
is little variation between bands. 

N(m) is the mean number of galaxies per magnitude 
per field area, which we derive in Section lSTTl The magnitude 
limits of the survey are given by mbright and mfain t . 

I n the case of the 1/V^ax luminosity estimator dSchmidti 
Il968h . the inverse of the total completeness, {TF{0i, mi))~^, 
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Figure 3. Actual survey completeness on the sky, R{0), for each of the KHJrpbj samples used in this paper. 
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is used to weight the inverse ma x imum accessible volume 
of eac h galaxy i. F or the STY|Sa ndagc. Tamman n fc Yahill 
I^197fl^ and SWML llEfstathiouTEUis fc PetersonlllQSSh estT 
mators, TpiOjin) is re-expressed in terms of absolute mag- 
nitude M for each galaxy using 



M = m - 51og[di(2:i)] - 25 - k{zi). 



(7) 



Here, dL{zi) is the luminosity distance in Mpc, k{zi) is the A;- 
correction and Zi the redshift of the ith galaxy. Equation (O 
is written in a form that ignores the evolutionary corrections 
that become significant at higher redshift. Tf(0,M) then 
weights the number density for all M. 



3 LUMINOSITY FUNCTION PRELIMINARIES 

3.1 Number Counts and Normalisation 

Large scale structure plays an important role in defining the 
overall normalisation of any LF, and is a crucial factor in 
determining the optimal depth and sky coverage of the sur- 
vey. Arguably it is the single biggest contributor to discrep- 
ancies between recent LFs and related quantities (see e.g. 
[Cole et al. 2001: N orbcrg ct al. 2002; Bell ct al. 2003; D riVeil 
I2OO4I: I Driver et al.ll2005fl . While recent luminosity function 
determinations have been able to exploit surveys of unprece- 
dented depth and coverage, cosmic struct ure is alway s an 
issue for the faintest galaxies of any survey. ICole et al] esti- 
mate that cosmic variance could introduce systematic errors 
as large as 15 percent in number counts from their sample 
combining 2MAS S NIR photometry with 2dFGRS redshifts. 
iBeUeiaLl (|2003) has found the SDSS Early Data Release to 
be 8 percent overdense with respect to the 2MASS extended 
source counts over the whole sky, and the 2dFGRS to be 
slightly under-dense. 

To date, the 6dFGS has covered about a quarter of the 
sky, or ~ 2 to 3 sr depending on passband. Compared to re- 
cent near-infrared (if-band) samples, th e 6dFGS cu r rently 
cov ers about ~ 15 to 20 times the area of lCole et al.l ll200lll 
andlBell et alj | 2003h and around 30 percent more that of 
iKochanek et al .| j200lD . In terms of sample size, the 6dFGS 
is ~ 10 — 1 5 tim es larger than Kochanok ct al. (2001) arid 
lBell_et al] (|200^ and a factor of 4 greater than ICole et alJ 
ll20^F ~At optical magnitudes, current 6dF GS sky cover- 
age q uadruples that used for t he 2d FGRS iNorberg et alJ 
f2002) and SDSS lBlanton et alJ \2Q0± luminosity functions, 
although sample size is roughly half. When complete, the sky 
coverage of 6dFGS will be almost double its present extent. 

Figure 2] shows the differential number counts for 
6dFGS sources after correction for Galactic extinction 
fSect. IX^ . The roll-over near the survey limit is the result 
of the extinction corrections blurring the sharp survey cut- 
offs imposed on the original uncorrected magnitudes. Unlike 
imaging number count surveys, we have been able to use our 
redshift information to remove stellar contamination from 
the counts. The counts were averaged over the effective area 
of the 6dFGS, obtained by scaling the eventual survey area 
by the fraction of targets in hand. Table (Ogives our linear 
least-squares fit parameters a, b for the form 



N{m) = 10' 



am+b 



(8) 



o 
I 



13.00 



b. „_= 16.75' 



-3 -2 
Total Magnitude 



-1 
Survey Limit 



where m is magnitude and N{m) is the differential surface 



Figure 4. Differential number counts for the KHJrFbj samples 
of the 6dFGS. The Euclidean slope of 0.6 has been subtracted for 
emphasis. Magnitudes have been corrected for Galactic extinction 
and are expressed as an offset from each survey limit (indicated). 
The solid lines show the straight lines of best fit given in Table ISl 
Poisson errors are shown. 



density ( mag~^ deg~^). The counts are flat in all bands 
with the exception of K, (and to a lesser extent, rp). This 
rise has also been see n in the 2MASS number counts pub- 
lished previously (e.g. ICole et alJl200ll : Ijarretdliooi) . and 
was initially thought to be due to the well-known upwards 
bias o f number count data in regimes of dec lining signal-to- 
noise jMurdoch. Crawford fc Jauncevlll973l) . We tested this 
by generating artiflcial random samples with the same size, 
flux and error distribution of the real A'-band sample. This 
Monte-Carlo approach revealed that even uncertainties as 
large as AK ~ 0.2 were unable to raise a nominal Euclidean 
0.6 slope by more than 1 percent. Typical 2MASS flux er- 
rors are closer to AK ~ 0.1 at K = 12.5. It is possible 
that the rise in K reflects the fact that the fainter bins are 
capturing more of the intrinsically lower surface brightness 
late-type galaxies, although no such rise is seen in the higher 
signal-to-noise H or J counts. 

Our source density of 3.2 galaxies mag~^ deg_^at_K_ = 
12.5 is com parable to the 3 m ag~^ deg~^ found bv lBell et alJ 
J2003h and lCole etaP J200ll) for the 2MASS XSC in Kron 
magnitudes. These are typically 0.1 to 0.2 m ag brighter 
than the original total mag nitudes of 2MASS llCole et alJ 
l200ll) . although it should be noted that the total mag- 
nitudes we now use for 6dFGS are the recently revised 
ones (Jarrett, priv. comm.). The 6dFGS number counts 
are identical to 2MASS number counts at Kiso = 12.2 
Jjarretdliooi) . allowing for the fact that 2MASS isopho- 
tal magnitudes are typ ically 0.1 to 0.2 mags fainter than 
total magnitudes. Like iBell et all we find that the Sloan 
Early Data Release counts arc somewhat higher than av- 
erage (3.5 mag~^ deg~^; IBcI I ct al. 2003), the likely re- 
sult of an overdensity in its comparatively smaller sur- 
vey region (370 deg'^). In the optical, the 6dFGS 6j-band 
counts of 2.9 mag~^ deg"'^ (6j — 16.0) are closer to 
the 2dFGRS cou nts in the S DSS -2dFGRS overlap region 
(2.9 mag~^deg~^: lNorberg'et al. 2002), although the NGP 
and SGP 2dFGRS regions (3.3 mag~^ deg~^) remain within 
the errors. Ultimately, such comparisons are limited by the 
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Table 3. Best fit parameters to extinction-corrected 6dFGS 
number count data. 



Band 


Fit 


range 


Number of 


Fit Parameters 




(mag) 


(mag) 


sources 


a 


b 


K 


10.75 


12.35 


113988 


0.636 


-7.132 


H 


10.95 


12.65 


90317 


0.588 


-6.769 


J 


11.75 


13.40 


93831 


0.583 


-7.114 


Tf 


13.60 


14.90 


64043 


0.621 


-8.782 




14.70 


16.10 


66833 


0.577 


-8.769 



systematics of different galaxy photometry techniques and 
magnitude definitions. 



3.2 Disentangling Hubble and Peculiar Motions 

The line-of-sight velocity measured for a galaxy has two 
components: the Hubble flow due to the expansion of the 
universe, and the peculiar velocity due to the combined grav- 
itational attraction of all large-scale structures. To infer true 
luminosity distances from the Hubble flow velocity, (assum- 
ing z < l: lHubbldlT93i) . we must flrst remove the pecu- 
liar motion component, especially given the low redshifts 
for much of our sample. 

We have adopted field fiow corrections from software 
written by J. P. Huchra f or the HST Key P roject to mea- 
sure Ho (Appendix A of iMould et~ai] l2000l) . This uses a 
simple linear fiow field with three spherical mass concen- 
trations representing Virgo {cz = 1035kms~^), the Great 
Attractor (GA; 4600 kms~^) and the Shapley Supercluster 
(SSC; 13800 kms~ ^). We also examined the fl ow models of 
iTonrv et all feOOd) and iBurstein et all (Il989l) . but did not 
use them since they neglect the SSC, which is a significant 
infiuence on much of the 6dFGS volume. 

Peculiar velocity corrections are most significant for 
cz < 10000 kms"\ Fi gure|K| shows how removal of the pe- 
culiar motion alters the absolute magnitudes for galaxies in 
the direction of the GA and SSC. The Virgo cluster direc- 
tion is not shown, as it lies beyond the northern limits of 
the 6dFGS. In the vicinity of the GA, galaxy luminosities 
are affected by as much as ±0.2 mags while for the SSC the 
effect is clearly less (< 0.1 mag). For galaxies within either 
mass concentration, the correction can be almost a magni- 
tude. This is because the fiow model sets the galaxy velocity 
to that of the mass concentration if the galaxy lies within 
its sphere of infiuence. 

With median redshifts in excess of 20 000kms-\ most 
recent surveys do no t include corrections for peculiar veloci- 
ties. An exception is|KpchMLek_et_al] l|20 0l[l, who employed 
the local fiow model of l^nn^^Lri|2000l) . lKochanek et alJ 
also employ a low-cz cut-off of 2000 kms~^, so as to exclude 
those local galaxies with a significant dependence on fiow 
corrections, at the expense of 5 percent of their sample and 
~ 1.5 mags at their faint-end limit. 

The nominal 6dFGS low-cz cut-off is 600kms~^, set 
so as to separate the obvious peak in Galactic sources at 
cz ~ Okms~^ from the rising redshift distribution of true 
extragalactic objects (Fig. (SJ. However, applying the flow 
model corrections broadens the distribution of all veloci- 
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Figure 5. Change in magnitude resulting from peculiar velocity 
corrections in two key sky directions (the Great Attractor and 
Shapley Supercluster, a and fe), and for the entire 6dFGS sam- 
ple (c). These show the difference between absolute magnitudes 
derived from the fiow-model-corrected cz values and uncorrected 
cz values. The diagonal bands are galaxies that the model has 
identified with the Great Attractor (4 600 km s~^) or Shapley Su- 
percluster {r^ 13 800 km s~^) mass concentrations, and ascribed 
the corresponding systemic velocity. 



ties, including those sources that are quite clearly Galactic. 
Therefore, we adjust the low-cz limit to 750kms~^ and use 
this throughout as our minimum acceptable redshift cut-off. 
Galaxies above the cut-off are used with their flow-corrected 
velocities from the model. 



3.3 Effect of Magnitude Errors 

The 6dFGS has been able to take advantage of the emerging 
availability of quality photometry over large areas of the sky. 
This is particularly the case for the NIR, where the digital 
2MASS photometry of the whole sky permits a uniformity 
and precision unprecedented in wide area redshift surveys. 
With this kind of precision, we are motivated to examine 
whether photometric errors modify the shape of the lumi- 
nosity distribution, and if so, by how much. For example, a 
non-gaussian tail in the error distribution can influence the 
bright end of the LF due to its rapid decline in this region. 

Figure |7| shows 2MASS photometric errors for JHK 
for 6dFGS targets. Typical errors are ~ 0.1 mag in K and 
~ 0.06 to 0.08 mag in J and H. In terms of absolute mag- 
nitude, these translate fairly uniformly across the range of 
observed luminosities (Fig. |HJ. Because of this, we apply 
the mean photometric error {a — 0.108, 0.083 and 0.065 for 
KHJ respectively) to all luminosities equally when com- 
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Figure 6. Effect of flow model on low cz values in relation to 
the low-c2 cut-off. The points show observed and model velocities 
relative to the lower and left-hand axes. The histograms show how 
the distribution projects along the same directions relative to the 
vertical scales shown on the right and above. The dashed lines 
indicate our \ow-cz cut-off for both observed (600kms~^) and 
model-corrected cz values (750kms~^). 




Magnitude 

Figure 7. Distribution of JHK photometry uncertainties for 
the revised 2MASS total magnitudes (points). The large solid 
blue points and long-dashed blue line both trace the la deviation 
of the uncertainties when grouped into 0.5-mag bins. The vertical 
solid red curves are half-gaussian representations of the same la 
deviation. The vertical short-dashed lines indicate the nominal 
6dFGS magnitude limits in each band. There are a small number 
of galaxies fainter than this limit due to the improvements to 
photometry after target selection. 



Figure 8. Distribution of JHK photometric errors as a func- 
tion of absolute total magnitude. The vertical lines are the result 
of dividing galaxies into 0.22-mag bins and co-adding gaussians 
equal to the la magnitude errors of the relevant galaxies. The 
short dashed lines indicate the ±lcr deviation and the long dashed 
histogram shows the distribution of 6dFGS sources on the same 
scale of absolute magnitude. The right-hand panels show the su- 
perposition of magnitude distributions for each bin (dotted lines) 
with the average for all shown by the solid line. The a for this 
distribution is indicated uppermost in each panel. 

puting its effect on the LF. Uncertainties in 6,, and rp are 
typically ~ 0.1, comparable with K. 

The functional form most commonly used to match the 
galaxy luminosity distribution is that proposed bv lSchechted 
il97(il . which gives the space density of galaxies per unit 
absolute magnitude, (j){M), as 

('1q0.4(M. -M)\a+1 

m) rfM = 0.4 In 10 0. ^^^p[,oo.4(M.-M)] (9) 

The parameter 0, is a normalisation constant, while a is 
the slope for faint galaxies and M* is a characteristic mag- 
nitude at the transition from the exponential to power law 
dependence. While both t he Schechter Function an d Press- 
Schechter Mass Function iPress fc SchechteJll974h are su- 
perficially similar, they are not related in any simple phys- 
ical way. The exponential cutoff at the bright end of the 
Schec hter Function is due to cooling (e.g. Rocs & Ostrike^ 
Il977^ and perhaps feedback by AGN. If it were caused only 
by the high mass cutoff of the Press-Schechter Function then 
Mt would be ~ 4 magnitudes brighter. 

In order to explore the effect of magnitude errors on this 
analytic fit, we convolved Schechter functions across a range 
of a values with the gaussian error kernels discussed above. 
Figure |3] shows how a generic Schechter function changes in 
log (^i-space when convolved with gaussian magnitude distri- 
butions of £7 = 0.07 to 0.1. As expected, only the brightest 
luminosities in the region M < M« are affected. Perhaps 
more surprising is that the amount of change is highly in- 
sensitive to the faint-end slope a for typical values. In the 
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Figure 9. The open circles sliow tlie offset between an ordinary 
Schechter function (log <f>) and one convolved with the gaussian 
error kernels for KHJrpbj (log t/iconv ) • The result is independent 
of a. The solid lines show the functional form of Eqn. llUi with 
V = 7.37 for Krpbj, 9.15 for H, and 12.25 for J. The vertical 
dotted line marks the location of M« . 

case of J, the magnitude errors are too small to have any 
appreciable effect on the Schechter fit. In the interests of 
computational efficiency we have used the empirical rela- 
tionship 

log (pconv = log + gA^.-M-u ^^q^ 

with values of ^; = 7.37 to match the convolution offsets in 
Krpbj, 9.15 for H, and 12.25 for J. Figure shows each of 
these curves. These offsets were then applied when finding 
optimal Schechter function fits in all bands (Section 2J. 

3.4 Extinction and Cosmological Corrections 

One of the primary science drivers for the NIR selec- 
tion of the main 6dFGS samples was the low impact of 
Galactic extinction over much of the sky . We use the 
maps of ISchleeel. Finkbeiner fc Davij il998l) to correct our 
total magnitudes using the relative extinction values of 
A/Av = 0.112,0.176,0.276,0.810,1.236 for KHJr^h, re- 
spectively. Typically this furnishes extinction corrections of 
up to ~ 0.04 mag in K, ~ 0.1 mag in J, and ~ 0.3 in 6j (re- 
call that the respective Galactic latitude limits are | & | > 10° 
for JHK and |6| > 20° for bar-p). Extinction corrections 
were applied after the initial apparent magnitude selection. 

Our cosmological fc-c orrections use the values com- 
puted by iPoggiantil |W97|), linearly interpolated with red- 
shift. Given the GdFGS selection bias towards early-types 
we adopt her fc-corrections for an elliptical galaxy with solar 
metallicity and exponentially decreasing star-formation rate 
with e-folding time of 1 Gyr. These corrections are within 
0.1 mag of those used bv iKoch ane^^^alJ (|20Mi) for t heir 
A'-band sample, and also those of lNorberge^dT(^2002^ for 
the 2dFGRS &j-band sample. 

We do not app ly evoluti onary e-corrections (e.g. 
iNorberg et~. 2002; Blanton et al. 2003'). since the bulk 
of our galaxies are seen at relatively recent cosmic time 
(z<0.1). 



4 LUMINOSITY FUNCTIONS 
4.1 6dF Galaxy Survey 

Several techniques have be en devised over the years for mea- 
suring galaxy LFs (see e.g. I WillmeJI 1 99'3) . We have adopted 
the three most commonly used — 1/Knax, STY and SWML 
— all of which have merits and sho rtcomings. The 1/Vmax 
method llSchmidtlll96l: lFeltenlll973) has the advantages of 
simplicity and no a priori assumption of a functional form 
for the luminosity distribution. Unlike the others, it also 
yields a fully normalised solution and is insensitive to mag- 
nitude incompleteness. Its main disadvantage is its assump- 
tion of a homogeneous, unclustered galaxy distribution. 
The STY method ( Sandagc, Tammann fc Yahil 1979) over- 
comes these difficulties, being a maximum-likelihood estima- 
tor with the added advantage that it does not require bin- 
ning. However, it presupposes some explicit functional model 
for the luminosity function. The Stop- Wise Maximum Like- 
hhood method (SWML; , Efstathiou, Elli s fc Peterson 198^) 
has similar advantages to the STY method, but does not 
assume an explicit form for the luminosity function. It pro- 
duces a binned maximum-likelihood estimate of the LF. 
STY and SWML provide shape information but not the 
overall normalisation in the manner of 1/Vmax. We refer the 
interested reader to the individual references for detailed 
methodologies. 

We have applied the 1/Vmax, STY and SWML meth- 
ods to each of the 6dFGS samples in KHJrpbj. The STY 
and SWML LFs were normalised by a minimisation with 
respect to the equivalent 1/Vmax distribution. 

Sources were excluded if they were outside the red- 
shift range 0.0025 < z < 0.2 or had total completeness 
TF{0i,mi) < 0.6. Completeness corrections had much less 
impact on the 1/Vnax results compared to oth er meth- 
ods, a point that has been noted by others (e.g. iKaldarel 
l2001fl . Those sources with extinction-corrected original mag- 
nitudes outside the range mbright to mfaint were also ex- 
cluded. Here, mfaint is simply the nominal survey limit 
(12.75, 13.00, 13.75, 15.60, 16.75) corrected for extinction in 
the direction of the source, while mbright = 8.75, 9.0, 9.75, 
13.0, and 14.0 for KHJrFbj respectively. Particular care was 
taken to exclude sources with rp < 13 or foj < 14 because 
of saturated rp or b,^ photometry in the original Southern 
Sky Survey plate material. The sizes of the samples actually 
used in the analysis are given in Tabled 

Figure 1101 shows the A"// Jr p&j LFs derived for the 
6dFGS via these methods. Table 0] summarises the best- 
fitting error-convolved Schechter functions, while Tables |S] 
and |S| list the non-parametric 1/Vmax and SWML distribu- 
tions. All three methods are in excellent agreement, insofar 
as the Schechter functions represent the data. Moving to 
bluer passbands we find a slight steepening of the faint-end 
slope, although not to the extent that has been reported pre- 
viously for dense environments. This is a consequence of the 
optical samples harbouring a larger population of late-type 
galaxies, which are typically of lower luminosity. 

The large sample size of the 6dFGS, even at this inter- 
mediate stage of the survey, has allowed us to beat down ran- 
dom uncertainties over a considerable range of luminosity. 
As a consequence, we can clearly see that a Schechter func- 
tion (Eqn.jnj is only an approximate match to the real galaxy 
luminosity distribution in all bands (upper panels in Fig- 
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Figure 10. Luminosity functions for the 6dF Galaxy Survey, derived from the 1/Vmax (green open circles), SWML (red solid circles) 
and STY methods (blue dashed curve). The inset shows the 1, 2 and 3fT confidence contours of the STY fit. The upper panel shows the 
1/Vmax and SWML residuals relative to STY (i.e. the deviations from the best-fitting Schechter function). 
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Table 4. Schechter function parameters for the STY luminosity functions of 6dFGS. The convolution parameter v used for Eqn. 1101 
is given for each band. 



Sample 


Fit Range 


M, - 5 log /i 


a 




V 




(mag) 


(mag) 




(h^ Mpc-3) 




K 


[-15.50,-28.85] 


-23.83 ± 0.03 


-1.16 ±0.04 


-2.126 ± 0.005 


7.37 


H 


[-16.00,-28.50] 


-23.54 ±0.04 


-1.11 ±0.04 


-2.141 ± 0.005 


9.15 


J 


[-15.00,-27.50] 


-22.85 ± 0.04 


-1.10 ±0.04 


-2.148 ± 0.005 


12.25 


rF 


[-14.00,-25.00] 


-20.98 ± 0.05 


-1.21 ±0.04 


-2.081 ± 0.006 


7.37 


h, 


[-12.50,-24.00] 


-19.91 ± 0.05 


-1.21 ±0.05 


-1.983 ± 0.006 


7.37 




(e) 



(LF orig) - (LF with one-component compi fit) 



-18 -20 -22 -24 -26 -28 

- 5 log,;, h 

Figure 11. Same X-band completeness measurements as shown 
in (a) - (d) of Fig. |21 (top panel), except fit with a single- 
component exponential, (e) The difference between the lumi- 
nosity function using these models and that using the double- 
exponentials of Fig. 1^ for K. The change is negligible in all bins 
except the very faintest. 



ure llUH . Specifically, it is unable to turn over sliarply enough 
around M«, and tends to slightly underestimate the true 
space density of Af, galaxies while at the same time slightly 
overestimating that of M ~ (A'/* - 1.5) and M > {M, + 1.5) 
systems. The differences are typically ~ 15 percent around 
Mt, and 30 to 40 percent beyond. Similar outc omes have 
been found in the LF fitt ing of deep SDSS data iBell et alJ 
12003 *: 'Blanton ct al. 2005") , prompting these authors to use 
double Schechter functions, or a hybrid Schechter function 
plus power law fit. 

In Sect. 12.31 we claimed that the luminosity function is 
unaffected by small deviations from the model used to fit 
the magnitude completeness. Figure [TTI shows the result of 
recomputing the if -band luminosity function using a single- 
compo nent exponential, such as those used bv lCoUess et alJ 
('2001') for 2dFGRS. We find there is negligible change in the 
luminosity function (Fig. Illfe ). which appears to be fairly 
insensitive to the exact form used to model magnitude com- 
pleteness. 
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Figure 12. (a) Apparent magnitude-redshift relation for the 
most luminous members of the if-band selected sample (Mjf < 
—26.5). (fe) Absolute magnitudes in KHJrwbj for the same Mk < 
—26.5 sample. Also shown are the Mt, values in each band (open 
circles and dashed line). Cosmological fc-corrections have been 
applied. 



4.2 Luminosity Function Extremities 

At faint luminosities M > (Af« + 2) we find marginal evi- 
dence for an upturn ^ as has been suggested by a number of 
authors for clusters (lDriveJ2004l . and references therein) . At 
the bright end we see a prominent upturn at A/ < (Af, — 2), 
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particularly in the K and fej samples, as has been noted by 
several authors previously. This excess of luminous objects is 
presumably due to brightest cluster galaxies, which are pro- 
duced by the special merger and accretion processes that 
come into effect in the high density regime at the centre of 
cluster gravitational potentials. We experi mented with fit- 
ting t his feature using the gaussian LF of ISaunders et all 
l|199Cl their Eqn. 6.1), but it only improved the bright-end 
fit at the expense of poorer fitting elsewhere. 

We studied the luminous galaxies responsible for the 
bright-end upturn in detail. Incorrect redshifts or the pres- 
ence of a non-gaussian tail on the magnitude error distri- 
bution could also cause an upturn if present in our data. 
Figure ^Ja) shows the magnitude-redshift relation for the 
most luminous galaxies with Mk < —26.5. It demonstrates 
that they span a range of of redshift and brightness, as ex- 
pected for a luminous sample that should be visible over 
a large survey volume. Figure I12f b') shows the fc-corrected 
absolute magnitudes in KHJrpb.^ for the same sample of 
galaxies. Some galaxies do not have matching magnitudes 
in the shorter wavelength bands due to non-selection by the 
original apparent magnitude cut. Also shown in Fig. ll2f fc') is 
the location of M, in each band. From this we observe that 
it is the same luminous galaxies that are responsible for the 
bright-end upturn in all bands. Hence, galaxy membership 
of the luminous bins is real and not the result of an extended 
tail in the magnitude error distribution. 

We also paid special attention to the spectra and red- 
shifts of this luminous sample. Figure lT^ shows a sub-sample 
of these with spectra and bjK imaging obtained from the 
6dFGS database. Generally, the spectra have signal-to-noise 
ratios that vary from excellent to poor. This is not sur- 
prising given the large range of apparent magnitude seen 
in Fig. \12i a). The imaging shows that there are a num- 
ber of galaxy pairs in this luminous sample. While they are 
not close enough to have inflated magnitudes due to im- 
age mergers — 2MASS can individually identify sources 5 " 
apart — it supports the idea that many of these galaxies 
inhabit dense envi ronments. In terms of redshift quality, Q 
J Jones et al.ll2004l) . around half the sample have Q = 4 in- 
dicating that the redshift measurements are secure. Of the 
remainder, half have Q — 3 while the other redshifts come 
from the literature. We can therefore rule out spurious red- 
shifts as a cause of the bright-end upturn. 

4.3 Comparison to Other Surveys 

Figure 1141 shows how the 6dFGS LFs compare to those of 
other recent surveys. In the K-hand, the 6dFGS LF agrees 
remarkab ly well with the earlier survey s using 2MASS pho- 
tometrv JCole et al ."200l': 'Koc hanek et alJl200ll: iBeU et alJ 
l2003l:lEke et alj2005il . .Eke et al.1 find a faint-end slope which 
is flatter than ot her recent LFs, including ours. The original 
iKochanek et al.l magnitudes have had 0.05 and 0.135 mags 
subtracted to convert from isophotal to Kron magni tudes, 
and then to our total magnitudes. The ICole et all Kron 
magnitudes have also had 0.135 mag subtracted to convert 
to total magnitudes. The larger sample and sky area of 
the 6dFGS provide ~ 1.5 to 2 mags additional coverage at 
both the bright and faint ends of the LF. We find that all 
previous surveys have tended to slightly underestimate the 
bright end, the likely consequence of smaller survey areas. 



iKochanek et alj find a steeper faint-end which they attribute 
to higher completeness and their brighter apparent magni- 
tude limit reducing sources of systematic error. However, 
none of the other K-h and surveys, including the 6dFGS, 
are as steep. iBell et all claim to see evidence of a faint-end 
upturn in their K-hand LF and use an ordinary power law to 
fit the faint end beyond Mk — 5 log h — —21. The faint end 
slope of our ordinary Schechter fit is a reasonable match to 
theirs over this range, and supports their estimated surface 
brightness incompleteness correction. 

The LFs i n Fig. ^^b) show how the lCole et al] ||200J) 
andtek^^y ( 2005) LFs compare with 6dFGS in J. Again, 
the ICole et aO Kron magnitudes have been transformed to 
total magnitudes by subtracting 0.135 mag. As with K, the 
J-band luminosity functions are in close agreement over 
those luminosities in common, although the 6dFGS yields 
slightly h igher densities of high-luminosity sources. The 
lEke et aP LF has a fiatter a and slightly higher normali- 
sation than 6dFGS. 

Figure I14f c') compares the 6dFGS rp-band LF with 
those from the Las Campanas Redshift Survey (LCRS; 
|LirLet_aD [l99 d) an d Century Survey with photographic 
JCeller et alj Il997ft and CCD photometry llBrown et alJ 
|20^F ! The 7?-band magnitudes of these surveys have 
been shifted to rp by subtrac t ing . 1 mag , following 
iFukugita. S himasa ku fc Ichikawal il995l) for local galax- 
ies of inter mediate type. We also show the deep LF of 
iBlanton et ali (2005) for a special low- luminosity subset of 
SDSS, which wc have transformed to rp from SDSS r by 
subtracting 0.34 mag. Our rp-band LF most closely agrees 
with the LCRS and SDSS surveys, and we suspect that 
cosmic variance is largely responsi ble for differen ces with 
the much smaller Century Survey. iBlanton et al.l see evi- 
dence for an upturn in the faint end and adopt a double- 
Sc hechter functio n to match it. Although not as pronounced 
as IBlanton et all the 6dFGS LF faint end shows marginal 
evidence of a steeper rise fainter than M — 51og/i —19, 
although this is by no means conclusive. Typically, a ris- 
ing faint end is associated with dense cluster environments 
jPriveil I2OO4I and references therein). The IBlanton et al] 
sample was deliberately chosen to be low luminosity, and 
therefore, of comparatively small volume. It is possible that 
the upturn they observe is the result of an overdensity in a 
nearby portion of their sample volume. 

In Fig. inid) we show how the 6 dFGS bj-band 
LFs compare to thos e from the 2dFGRS jNorberg et alJ 
l2002h . Stromlo- APM iLovedav et alJll992ll and ESO Slice 
surveys (ESP; IZucca et^^Tl^gTl) . We have also trans- 
formed the SDSS .g-band LFs of iBeU et alJ J2003l) and 
^lanton ct al. (2005) by adding 0.25 mag to match 6j, fol- 
lowing IBlanton et all We also sh ow the LF from th e Mil- 
lennium Galaxy Catalogue (MGC: lbriver et al.l2005ll which 
has been transformed from their bl ue passband to ours 
through bj = Bmgc — 0.13. The LF of IZucca et al.l has been 
modified to take into account the cosmological model we 
have adopted. 

Figure ITU shows that the 6dFGS faint end slope {a = 
-1.21) is closest to that of 2dFGRS and ESP, although 
the 6dFGS normalisation is lower than the others, possi- 
bly due to our non-use of evolution corrections (which are 
most significa nt in bluer pas s bands). Th ese magnitude cor- 
rections (e.g. iNorberg et all iBell et all) tend to lower the 
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Figure 13. Spectra of some luminous upturn sources at Mk < —26.5. Tlie labeling of spectral features indicates position only and 
does not imply detection of the feature. Also shown are postage stamp images in bj and K, taken from SuperCOSMOS and 2MASS 
respectively. North is up, east is left, and image size is indicated in the lower right-hand corner. 
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Table 5. Near-infrared l/Vmax and SWML luminosity functions from the 6dFGS. 
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Figure 14. Comparison between tlie 6dFGS luminosity functions for KJr^bj and tliose of other surveys. Only the STY and SWML 
fits to 6dFGS have been reproduced. Ordinary Schechter functions are shown throughout except for the double Schechter function of 
Blantoi^talJ (c and d) and the composite Schechter-and-power-law form of iBell et alJ for K (a). The Ka.nd J Kro n magnitudes of 
Col^^^aLf have been converted to the total magnitudes used here by subtracting 0.135 mag; the isopho tal magnitudes of lKochanek et al] 
hav e similarly been corrected by subtract ing 0.185 mag. The passband transformations rp = r — 0.34 ffor lBlanton et all) an d r-p = f?— 0.10 
(forlLii^^a^ Geller et alilBrown et all) have been used in (c). In (d), bj = q + 0.25 has been used to transform the LFs of lBlanton et all 
and iBell et al. while fej = Bmgc ~ 0.13 was used on the LF of iDriver et aI.L The approximate absolute magnitude coverage of each 
survey is indicated by the extent of its curve. 



bright end end relative to the faint end, and lBlanton et all 
has demonstrated ~ 20 to 30 percent variations are possi- 
ble depending on the strength of the applied correction. 
We note that the 6dFGS norma lisation (log(^, = —1.983) 
most closely matches that of the lBlanton et all SDSS sam- 
ple {log(j)t = —1.928), which (like the 6dFGS) was not cor- 
rected for evolution. We do not find evi dence for the steep 
faint end upturn found bv lBlanton et all despite having suf- 
ficient luminosity coverage to test for it. We suspect a local 
overdensity is the likely culprit for this feature of their data. 
Ultimately, subtle differences are inevitable because of the 
way deeper but narrower surveys (SDSS, 2dFGRS, ESP) 



sample the very nearest large-scale structures compared to 
the shallower but wider surveys such as Stromlo-APM and 
the 6dFGS. 



5 LUMINOSITY DENSITY 

The total luminosity per unit volume, or luminosity density, 
j, is a definitive observable of any galaxy population. It is the 
integral of the luminosity function, weighted by luminosity, 
and is given in solar luminosity units by 

j = 0, ■ io-O-4(A^--A'^0) . r(a + 2). (11) 
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Figure 15. Integrated luminosity density for 6dFGS across 
bjr^JHK in units of 10* Lq h Mpc~^ (large red solid circles, solid 
red line) . Comparitive v alues from other recent surveys are also 
shown. Both K values of lBell et al.l l2003^ are included: one with 
and one without correction for low surface brightness selection 
ef fects. We have also reproduced the spectral energy distribution 
of iBell et a l. ( 2003) for a 12Gyr old stellar population with ex- 
ponentially decreasing star formation rate, with e-folding time of 
4Gyr (small solid circles, dashed line). 



Figure 1151 shows the luminosity densities calculated from 
Eqn. CIJ and the 6dFGS STY fit parameters in Ta- 
ble 3] We have used solar absolute magnitudes of Mq = 
(5.442,4.447,3.660,3.319,3.280) for h.r^JHK respectively, 
made avai lable by C . Willm er^ . We note that in the simi- 
lar plot of iBell et alJ (|200d) (their Fig. 15) the &j point for 
iNorberg et alJ l^2002^ is located ^ 700 A blueward of its true 
value. 

Figure ITKl also shows how the luminosity densities from 
GdFGS compare to those of other recent sur veys. Although 
much has been made (see e.g. IWright 200ll; iNorberg et alJ 
|2002|) . of ~ 2 x discrepancies betw een the early SDSS opti- 
cal estimates llBlanton et alJl200 J) and 2dFG RS in the NIR 
llCole et al.ll200lh . the revised SDSS values of iBlanton et alJ 
i2003r . corrected for cosmological evolution, have largely 
ameliorated these differences. The 6dFGS values agree with 
most recent estimates, finding a K-ha,nd value at the lower 
end of recent results for this band. This is also the expecta- 
tion of the spectral energy distribution (SED) for a 12 Gyr 
old stellar population with a 4 Gyr e-folding exponentially 
decre asing star formation rate, (dashed line; after Bell et al. 
|20o3). Models with more rapidly declining star formation 
rates (r^ 2 Gyr) can produce SEDs more luminous in K, 
but at the expense of overestimating optical luminosities by 
factors of 2 or more. 

While the luminosity density is a convenient quantity 
through which to compare different surveys and infer a 



http: / /www. ucolick.org/ r^cnaw/sun.html 



volume-averaged SED, it ultimately relies on a fit to the 
LF and its extrapolation beyond accessible luminosities. We 
find that the STY fits consistently overestimate the inte- 
grated luminosity of the S WML method by ~ 1 to 3 percent 
across the observable range. The observed luminosity dis- 
tribution typically accounts for around 95 percent of total 
luminosity, with the shortfall due almost entirely to faint- 
end extrapolation. This is comparable in size to the extra 2 
to 5 percent integrated luminosity measured by the 6dFGS 
in reaching ~ 1 to 2 mags fainter in the NIR passbands. 



6 SUMMARY 

We have presented luminosity functions in KHJrpbj from 
the 6dF Galaxy Survey and derived luminosity densi- 
ties from them. The 6dFGS target samples are drawn 
from the 2MASS XSC and SuperCOSMOS catalogues 
and are near-complete to limits of (K , H , J ,rF,b,,) = 
(12.75,13.00,13.75,15.60,16.75). The samples used here 
represent around half the final 6dFGS in terms of sky cov- 
erage and sample size. The 6dFGS is already the largest 
NIR-selected redshift survey by more than an order of mag- 
nitude in both coverage and sample size. Compared to the 
SDSS and 2dFGRS optical LFs it has around four times the 
sky coverage and half the sample size. 

The 6dFGS luminosity functions have been calculated 
using the 1/Knax, STY and SWML estimators. The effects of 
magnitude and field-related incompleteness have been char- 
acterised and taken into account. Luminosity distances have 
been corrected for the effects of peculiar velocity, which can 
amount to differences as large as ~ ±0.5 mag for some super- 
cluster members. The influence of photometric errors on the 
shape of the LF has been taken into account for KHr^bj, 
but was of little consequence for J. 

These new 6dFGS LFs probe ~ 1 to 2 mag fainter in 
absolute magnitude than previous surveys in the NIR pass- 
bands and to comparable limits in rp and 6j. We obtain 
excellent agreement between the 1/Vmax and SWML esti- 
mates of the LFs, but find that a Schechter function (fitted 
by the optimal STY method) either under- or over-estimates 
these values by as much as 15 to 40 percent. While the for- 
mal uncertainties on our STY fits are typically a few percent 
or less in log(/)«, M* and a, we show that, in comparison to 
the observed LF, a Schechter function is unable to decline 
rapidly enough at the bright-end and remain as flat at the 
faint end. 

The 6dFGS results are generally in excellent agreement 
with other recent LF determinatio ns in the NIR. We do no t 
find as steep a faint-end slope as iKochanek et al.| i200ll) . 
and suspect this is due to their shallower depth and sub- 
sequent brighter faint-end limit. The 6dFGS rp-band LF 
most closely matches those of the Las Campanas Redshift 
Survey ^in ct al. 1996) and SDSS (Blanton et al. 2005), al- 
though we find o nly marginal ev idence for the faint-end up- 
turn claimed by iBlanton et alj In 6j, the 6dFGS LF has 
a nearly identical faint end slope to those obtained by 
2dFGRS f Norberg et al."2002') and the ESO Slice Project 
iZucca et a l. 1997), alt hough the 6dF GS normalisation is 
closer to that found bv IBlanton et alT Neither this survey 
nor 6dFGS used evolutionary corrections. Furthermore, we 
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see no evidence for a faint-end up-turn in any of the 6dFGS 
LFs. 

The luminosity densities derived in all five 6dFGS pass- 
bands concur with most other recent measurements, and 
support a if -band value at the lower range of recent values. 
This is consistent with the effective mean galaxy spectral 
energy distribution being represented by an old stellar pop- 
ulation with moderately decreasing star formation rate. 
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